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We formulate particle creation phenomena in Bose-Einstein condensates in terms of conserving 
gapless mean field theory for weakly interacting Bose gases. The particle creation spectrum is 
calculated by rediagonalizing the Bogoliubov-de Gennes (BdG) Hamiltonian in mean field theory. 
The conservation implies that quasiparticle creation is accompanied by quantum backreaction to the 
condensates. Particle creation in this mean field theory is found to be equivalent to that in quantum 
field theory (QFT) in curved spacetime. An expression is obtained for an effective metric affected 
by quantum backreaction. The formula for the particle creation spectrum obtained in terms of QFT 
in curved spacetime is shown to be the same as that given by rediagonalizing the BdG Hamiltonian. 
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I. INTRODUCTION 

Particle creation is one of the most important quantum 
phenomena predicted by quantum field theory (QFT) in 
curved spacetime; in particular, it has important implica- 
tions for cosmology and gravitational physics [l[. When 
spacetime evolves dynamically, matter field quanta in 
spacetime are created even when no quanta initially exist. 
Hawking radiation [2[ is the most well-known example of 
particle creation. When a black hole is formed by grav- 
itational collapse of a star, Hawking quanta are created 
and then emitted from the black hole as blackbody radi- 
ation. Hawking radiation is commonly believed to cause 
a backreaction to the black-hole spacetime, resulting in 
black-hole evaporation. However, a precise description of 
the backreaction of quantum particle creation to classical 
spacetime dynamics has not been developed yet |3[ (for 
recent progress see, for example, [J, [a]). 

A scalar-field wave propagating in curved spacetime is 
analogous to a sound wave propagating in a fluid [6[ . In 
this analogy, background fluid flow corresponds to curved 
spacetime and provides the spacetime metric. The phase 
fluctuation of a sound wave can then be considered a 
field in curved spacetime. This analogy enables us to use 
fluid systems to studymany field theoretical phenomena 
in curved spacetime 0, B ■ It is also possible to use the 
analogy to derive some implications of quantum backre- 
action to spacetime dynamics by investigating condensed 
matter systems [9l-tF3|. 

One advantage of using such an analogy is the possibil- 
ity of observing quantum effects experimentally. It is im- 
portant to consider quantum fluids to investigate quan- 
tum effects precisely. Bose-Einstein condensates (BECs) 
in trapped cold atoms |13l. [l4| are some of the best sys- 
tems in this context [15[. Acoustic black holes has been 
recently reported for the first time in BEC [ljj. One of 
the simplest methods for realizing a sonic horizon is just 
to expand a BEC [131 • The analogy of a BEC system 
is fascinating from a theoretical point of view because it 



provides a microscopic quantum description of quantum 
excitations. For example, a microscopic understanding 
of analog Hawking radiation has been obtained by inves- 
tigating it in BECs in Bogoliubov theory |18l - t2l| . 

Particle creation and its quantum backreaction to con- 
densates will occur in BEC systems and should be de- 
scribed by the theory of quantum fluctuations in con- 
densed matter systems, independently of any analogy 
to curved spacetime. Our aim is to formulate particle 
creation in condensed matter theory and to clarify the 
backreaction effect from quantum Bogoliubov quasipar- 
ticle creation to background condensates as well as its 
correction to particle creation itself. This is important 
because backreaction effects are expected to be observed 
in experiments because of particle number conservation. 

This aim will also be valuable for theoretical physics, 
if it succeeds in obtaining a method for dealing with the 
backreaction problem. For example, Schutzhold et al.[9| 
investigated quantum backreaction in dilute gas BECs, 
and pointed out that the effective-action technique does 
not yield the correct result in general. This indication 
is important for gravitational physics, because the pre- 
scription of quantum backreaction via the expectation 
value of the pseudo-energy momentum tensor T M „ to the 
Einstein equations, which is usually considered, in the 
form 



K'llV rydP-^T^ \^IIV/1 



(1) 



would not be sufficient. Therefore, the analogy provides a 
useful theoretical tool for considering the quantum aspect 
of the backreaction problem in QFT in curved spacetime. 
The quantum backreaction to condensates has been 
discussed apart from the analogy context. Braaten and 



Nieto [22|, and Andersen and Braaten |23( discussed 



semiclassical corrections to a general time-independent 
condensate using Hartree-Fock-Bogoliubov (HFB) the- 
ory, which is a particle-number-conserving mean field 
theory. Buljan et al. |2J] considered the backreaction 
of thermal quantum fluctuations to condensate dynam- 



ics using HFB theory. They performed a numerical 
simulation as an example of a backreaction in a finite- 
temperature BEC system. 

We consider backreaction effects of particle creation to 
condensates which are not included in these works |22| - 
l24j . Particle creation is expected to occur in systems 
with dynamical background field (condensates or curved 
spacetime, and so on). Furthermore, particle creation 
occurs because of a change in definition of particles [l[ , 
and can not be evaluated only by solving time dependent 
Gross-Pitacvskii (GP) and Bogoliubov-de Gennes (BdG) 
equations, even when these equations include backreac- 
tion terms as quasiparticle pair correlation. In the Bo- 
goliubov theory, particle creation is formulated in terms 
of rediagonalization of the BdG Hamiltonian [25|, l26f . 
representing change in definition of quasiparticles. Of 
course, in the Bogoliubov theory, the particle number of 
a condensate is conserved, and the backreaction to the 
condensate cannot be taken into account. 

HFB theory satisfies the conservation law within the 
total system. Therefore, HFB theory is a candidate the- 
ory for treating the backreaction problem. However, we 
do not consider this theory for the following reason: The 
backreaction of particle creation has two aspects. One 
is the backreaction to the condensates, and the other is 
the backreaction to the particle creation spectrum itself, 
caused by the correction of the condensates. We explore a 
consistent formulation including these two aspects. For 
the purpose of evaluating the second aspect, a gaplcss 
theory is necessary. In a gaplcss theory, by rediagonal- 
izing the BdG Hamiltonian, the condensate state can be 
considered in the calculation, and the mixing of the con- 
densate state and the redefined quasiparticle states can 
be evaluated. However, HFB theory is not gapless. Us- 
ing HFB theory, we can not evaluate backreaction cor- 
rections to the particle creation spectrum. Therefore, we 
do not consider HFB theory in the present article. There 
is another reason why HFB theory is not available. HFB 
theory has an energy gap in the excitation spectrum; this 
contradicts the Hugenholtz-Pines theorem [27| , which re- 
quires gapless excitation in the long- wavelength limit. 

Recently, Kita [28l - l3l| formulated an alternative mean 
field theory for BECs that both satisfies conservation 
laws and has gapless excitations. It is thus a promising 
theory for describing BEC systems. In this article, we 
consider the backreaction problem in terms of Kita the- 
ory. By use of this theory, we formulate particle creation 
that includes backreaction effects in BECs at zero tem- 
perature (i.e., we do not consider thermal fluctuations). 
The particle creation spectrum corrected by backreaction 
effects is obtained. We also formulate an analogy with 
curved spacetime in this theory. By considering the anal- 
ogy, we seek to obtain a model representing the quantum 
backreaction effect in the effective spacetime metric and 
to obtain some implications for QFT in curved spacetime. 

The article is organized as follows. In Sec. [Til we re- 
view some basic equations and eigenstates of the BdG 
Hamiltonian in Kita theory, which are necessary for the 



subsequent discussion. In Sec. IIII1 wc discuss particle 
creation based on QFT in curved spacetime. For that 
purpose, we construct QFT in a spacetime analog based 
on Kita theory. In subsection IIII A[ the correspondence 
between the equation of motion for a quasiparticle field 
and a quantum field in curved spacetime is clarified. In 
subsection IIII B[ a formula for the particle creation spec- 
trum is obtained using QFT in a spacetime analog. In 
Sec. IIV1 we formulate particle creation in BECs within 
BEC theory as rediagonalization of the BdG Hamilto- 
nian. The formula obtained for the particle creation 
spectrum is compared with that calculated in Sec. IIIII 
Section [V] gives the conclusions of the study. 



II. REVIEW OF BASIC EQUATIONS: BDG 
EIGENSTATES 

We consider weakly interacting bosons trapped by an 
external confining potential 14 x t(x). The interaction is 
assumed to be a contact interaction with a strength de- 
noted by Uq. The field operator ip can be decomposed 
into the condensate wave function ^ and the quasiparti- 
cle field 6 as 



i> = * 



(2) 



In time-dependent Kita theory, the equation for the con- 
densate wave function ^ is given by 



ih 



~dt 



L + C/ |*| 2 ) tt + 2C/ O (0 t 0)* + £/o( #)**, (3) 



where L = K + V cxt - /a. K = -|^;V 2 is the ki- 
netic energy operator and /x is the chemical potential. 
The extended Gross-Pitaevskii (GP) equation (Eq. (J3j)) 
includes quasiparticle pair correlations (<ft4>) and (4>(f>). 
These correlations are obtained from the renormalizcd 
Green's function [3l| and we assume that these quanti- 
ties are finite. The condensate wave function ty can be 
rewritten in terms of the condensate number density no 
and the phase S as: 



* 



n e 



is 



(4) 



The total number of condensate particles is given by 



N 



\V\ 2 d 3 x 



nod x. 



(5) 



In general, No is time dependent in this theory. 

The time-dependent BdG equation can be represented 
as 



ihd t & = H* 



where the BdG Hamiltonian is given by 



H 



Uo^ 2 - U, 



W + 2U 
-Uo^* 2 + U (tf tf) -W-2U 



(6) 



, (?) 



where we have defined W := L + 2it,oUq. $ is a two- 
component vector composed of 4> and (ft 



*(i,x) 



( 0(*,x) 
\<ft{t,*) 



(8) 



We introduce a similar vector for the condensate wave 
function, 



*(*,*) 



tt(i,x) 

vT(i,x) 



(9) 



It is found that 73 \I/ with the third Pauli matrix T3 is a 
solution of the BdG equation (Eq. (|6])), which implies 
that the excitation spectrum is gapless. 

Next, we review the eigenstates of the BdG Hamilto- 
nian for some quasi-static situations. We denote the 
eigenvectors of the BdG Hamiltonian as (f m , which have 
two components. They satisfy 



H<p m =E m (p m , 



(10) 



where to is an integer. As discussed in Ref. [29j, the 
orthonormal relation between eigenstates with positive 
eigenvalues can be assumed to be 



vVml^lVm') = $n 



(11) 



It can be shown that a positive eigenvalue Ej of the BdG 
equation with the eigenvector 



u, 



(12) 



is always accompanied by the negative eigenvalues — Ej 
with the eigenvector 

(13) 

Then, the normalization Eq. (|ll[) implies that 

( Uj \u k ) - {v k \vj) = S jk . (14) 

Equivalently, the eigenstates with negative eigenvalues 
satisfy 



{Vj\T3\V k ) 



Ojk- 



(15) 



In equilibrium, the BdG Hamiltonian has a zero eigen- 
value whose eigenvector is proportional to T3|^E') with \*f?) 
given by Eq. ©. We normalize this eigenstate as follows: 



¥0 '■= 



1 



=f 3 \I> = 



fo 



where we have defined the function tpo as 

1 



<Po := 



-A'. 



(16) 



(17) 



Noting that (*|*) = 2N , we find (volVo) = 1- 1* is 
orthogonal to the states |u,) as 



( u #3|<Po} =0= 



(18) 



As discussed in Ref. |29|, by noting Eq. ([14")) , we may 
assume 



giving 



(Uj\ip ) = {<Po\vj) = 0, 



(UjlVo) =°- 



(19) 



(20) 



For BECs in equilibrium, these eigenstates represent Bo- 
goliubov quasiparticles. 

Using these eigenfunctions, the quasiparticle field op- 
erator can be expanded as 



E 



Ik 



Uk 

-vt 



ll 



-Vk 



+(7o-7^)(_^). (2D 

The operators 7 m and 7^ (to, n = 0, 1, 2, • • • ) satisfy the 
boson commutation relation [7m7 7i] = &mn- The oper- 
ators 7fc and 7^ can be interpreted as annihilation and 
creation operators of quasiparticles, respectively. 70 can 
be interpreted as operators that annihilate a particle from 
the condensate, and 7g is the reverse operator. We have 

determined the minus sign before 7q in the last term of 
the expansion (Eq. (f2~Tj) ') in accordance with the discus- 
sion in Appendix |A1 



III. ANALOGY WITH QFT IN CURVED 
SPACETIME 



In this section, we discuss particle creation in BECs us- 
ing the analogy with QFT in curved spacetime. For this 
purpose, we construct QFT in curved spacetime within 
the conserving gapless mean field theory |28l - [3l| . First, 
we show that the equation of motion for quasiparticles 
can be viewed as an equation for a quantum field in 
curved spacetime. The metric of the spacetime analog 
includes the backreaction effect. Next, we obtain the 
particle creation formula based on QFT in curved space- 
time. 



A. field equation 

The quasiparticle field can be U(l) gauge transformed 
as 



:=e-' lS (h, (ft:=e iS (ft. 



(22) 



The extended GP equation (|3]) can be rewritten as the 
following two real equations: 



d t n + V • (n v ) - 2fl»a ((#) - (0t^t)) = 0, 



(23) 



hd t S + n -ft^o + V ex t — n + no?7o + jmu. 



1 2 



+2U (^4>) + \U (<#) + (0t0t)) = o, (24) 

where vo := KS7 S/m is the velocity of the condensate. 
The equations for and <^ are 

ihd t 4> = (L + 2n U + hd t S)4> + noUoft 

+2U {ftM-U {M>)&, (25) 

-ihdtft = (L* + 2n U + hd t S)ft + n U 4> 



+2U {ft&ft-Uo(ftft)j>, 


(26) 


where L is defined as 




L := K + V ext - u, 


(27) 


with 




K:=-^-(W + l (WS)f. 


(28) 



Performing a short calculation, we obtain 

K = K - ihv ■ V - -ih{V • vo) + -zmvl. (29) 

Density fluctuations and current density fluctuations are 
described by 

p' ■■= s/n ~(<f> + 4>' t ), j':=/5'vo + nov', (30) 

where we have introduced v' = V<&' with 

h 



$' 



2miJriQ 



(0 ~ f )- 



Taking the difference of Eqs. (|2~5j) and 
the following conservation law 

5 t/ 5' + V-j' = 0. 



(31) 
we obtain 

(32) 



In the hydrodynamic approximation, in which the den- 
sity gradient is smooth over the local healing length 
£ = h/y / 2mUono, the sum of Eqs. (|25|) and ([26)) gives 



UnX 



(d t + v • V) $', 



where the quantity X has been defined as 



X := 1 



2n n 



(33) 



(34) 



Note that X is real. Then, Eq. (|32|) can be rewritten as 



^g^dA' 



= 0. 



(35) 



In this equation, we have defined the symmetric tensor 
9iw as 



9V 



mc s ( -{cl-vl) ~v 0a 

UnX I -«06 Sab 



(36) 



where a, b = 1, 2, 3. The sound velocity can be described 
in terms of X as 



n U X 



(37) 



If we identify g^ v with a spacetime metric through ds 2 = 
g^dx^dx 1 ' , the field $' satisfies the equation of motion 
for a minimally coupled scalar field in curved spacetime 
with the metric given by Eq. (ptrJ)). It can thus be con- 
sidered as a field in curved spacetime. 

The effective metric includes contributions from quasi- 
particle pair correlations. If we drop the contributions 
and set X = 1, then the effective metric becomes the 
well-known metric without the backreaction (see, for ex- 
ample, [8j). We find that the correction to the effective 
metric appears in the form of UnX . However, the back- 
reaction effect is not included only in X: the quantities 
no and vn are affected by the backreaction effect through 
the pair correlation terms in the extended GP equations, 
Eqs. (|23 ]) and ([24 ]) . 



B. Particle creation spectrum 

We calculate the particle creation spectrum by using 
QFT in curved spacetime in the usual way, namely with- 
out considering a specific prescription for treating the 
backreaction to the background field. 

First, we introduce U(l) gauged wave functions Uj and 
ij, which satisfy the U(l) gauged BdG equation 



ihdt 



H 



(38) 



with the gauged Hamiltonian 



~ = W + 2Un{^4>) nnUn-Un{^ 

■ \ -mUn + Uoiftft) -W* - 2U (^4>) J ' [ ' 

where W := L + 2Uqiiq + hdtS. The normalized and 
gauged condensate wave function is defined by 



1 



Vo 



V2lVo 



"o- 



(40) 



Similarly, we define the gauged condensate wave function 
as ^ := e _ \& = y/rtQ. The normalization factor Nn in 
Eq. ([4TJ]) is the number of condensate particles at some 
time and it is a constant. It is found that 



T:>, 



Vo 



Vo 
-0o 



(41) 



also satisfies the gauged BdG equation with the Hamil- 
tonian Eq. ([55]). 

With these fields, the U(l) gauged field (f> — e~~ lS <f) can 
be expanded as 



5Z (w - l]vj) + 7o<Po - 7o^o 



(42) 



The scalar field <&' can be expanded in terms of wave 
functions 



fi("j + «j ) 



for positive j and 



/o = 



ImiJriQ 



(fio + tPa) 



mi^/21%' 



as 



^Eiu+^i 



(43) 



(44) 



(45) 



n=0 



The zero mode /o is a constant. 

Under the hydrodynamic approximation, it can be 
shown that the mode functions {/„} obey the equation 
of motion in the spacetime analog 



d u 



-99^dJ n 



0. 



(46) 



The details of this derivation are given in Appendix [Bj 

Next, we investigate the Klein- Gordon (KG) inner 
product between the mode functions. The KG product 
for functions F and G is defined as 



(F,G) 



KG 



I f {Fn^d^G* - GTn^F) Vdet hd 3 



(47) 



where n = fif'd^ is a unit normal vector orthogonal to 
the spacelike surface E; that is, n is a future-directed 
timclike vector. Additionally, Vdct hd 3 x is an invariant 
volume clement of the surface E, and h a b is the induced 
metric on E. If we choose a ^-constant surface as E, the 
induced metric is 



ds s = h a bdx a dx 



UoVX 



6 ab dx a dx 



(48) 



Then, the determinant of the induced metric is Vdet h 
(mno/UoX) 3 / 4 . The normal vector to the surface is 



i Uo 



1/4 



(dt + Vo-V). 



(49) 



For the constant /o, it is trivial to show that 

(/o,/oW=0. (50) 



Furthermore, 

(fj,fo)KG 



no 

2N 



(uj — v*)d 3 x 



(ip \Uj) - (Vj\<p ) 



(51) 



When the system is quasi-static, Eq. (fT9"]) can be applied 
to the preceding equation, giving 



C/j./o) 



KG 



0. 



(52) 



Therefore, KG products including the zero-mode func- 
tion /o vanish in quasi-static situations. The mode func- 
tions fj with positive j satisfy 



rfdj, = - 



n^UoX 



1/4 U X^ 

— (Uj-Vj). (53 

rn J 



See Appendix [B] for the detailed calculation. Using the 
orthonormal relation of Eq. (fl"4|) . it can be shown that the 
functions are normal each other, 



(fjtfk) 



KG 



-^ 



(54) 



Similar calculations show that {fT,f%)KG = ~m^J k an< ^ 
{fj,fk)KG = 0. Therefore, mode functions with positive 
j are normal each other. At the same time, they are not 
orthonormal. The mode functions should be orthonor- 
mal to identify the number of quanta in effective curved 
spacetime and the number of Bogoliubov quasiparticles. 
For this purpose, we normalize the field <&' as 



2iy?mno y H 

With the normalized mode functions 

Vh(uj +vj) 



/,' == 



li^JmriQ 



and 



.fo 



(tp + ip ) 



2iy/mno yT i^/mNo 

the normalized field tp can be expanded as 



v 



n=0 



(%fn + fnfn) 



(55) 



(56) 



(57) 



(58) 



The number operator of the field tp is given by 

N 3 =l]%, (59) 

which is the same as the number operator of Bogoliubov 
quasiparticles. We have exactly related one Bogoliubov 
quasiparticle with one quantum in the spacetime analog. 



As a simple example of particle creation, we consider a 
condensate that is quasi-static and has no quasiparticle 
excitation at initial time t = t\. For t\ < t < £2, the 
condensate evolves dynamically and becomes quasi-static 
again at t = t%. At the initial and final times t = ti 
(i=l,2), we expand the field as 



f 



OO 

E 



V In J n > In J n h 



(60) 



r(t) 



where {/A } is a complete set at each time t = ti. Each 



F (i) 



/„ consists of eigenfunctions of the BdG equations u,- , 



(0 



,(') 



,(*) 



w , and ip at each time t = ti. With annihilation and 

creation operators 7)1 and 7™ , we define one particle 
state of the quantum. 

The calculation of the particle creation spectrum based 
on QFT in effective spacetime is exactly the same as that 
given in [25| , except that the current expansion includes 
the zero mode /q . However, as in Eqs. (jSTJj) and ([52"]) . 
the zero mode does not contribute to the KG product, 
and does not alter the particle creation spectrum. There- 
fore, the particle creation spectrum is the same as that 
given in Ref. |25| ■ Explicitly, the expectation value of the 

a (n\ A (2)+ A (2) 

final number operator in the jth mode iv • := 7 7 



(i 



<0|JVf|0) (1) 



= E 



B 



./'■■ 



(61) 



where the coefficient Bjk is calculated as the KG product 
at the final time t = t2, 



R _ Cf(2) f (l) 



J KG 



d 3 *^ 2 ^* 



(2)* (l)s 



(62) 



This expression agrees with the one obtained in Ref. [25j , 
which did not consider any backreaction effect [38[ • How- 
ever, we should note that wave functions Uj, Vj, and tpo 
in Eq. (|62|) are solutions of the equations that already 
include backreaction effects. Therefore, although the ex- 
pression for the particle creation is the same, it differs 
quantitatively from the result obtained in Ref. |25j until 
the large condensate limit (Nq — > 00) is taken. 



IV. PARTICLE CREATION AS BDG 
REDIAGONALIZATION 

In this section, we consider particle creation within 
Kita theory as the rediagonalization of the BdG Hamil- 
tonian, without referring to any effective spacetime. 

Now, we consider the eigenspace of the BdG Hamilto- 
nian. The unit operator that acts on the eigenspace is 
assumed to be expressed in terms of eigenstates as 



E 
3 



l u i)( u #3- \Vj){Vj\h 



IVoXVol- 



(63) 



This complete set differs from that introduced in Ref. 
J29|, which includes a state that does not belong to the 
BdG eigenstates. Since we are interested in solutions to 
the BdG equation, we assume the complete set given by 
Eq. (|3]). 

If the condensate and quasiparticle system are com- 
pletely static, then no particle creation will occur. How- 
ever, particles will be created in general when the con- 
densate moves dynamically. As a simple example of par- 
ticle creation, we consider the same situation for BECs 
as considered in Sec. IIIIB1 namely the time evolution 
of condensate for t\ < t < t2 accompanied by the ini- 
tial and final quasi-static condensates. Then, the com- 
plete sets that diagonalizc the BdG Hamiltonian at t = t\ 
and t = t2 will be different in general, and annihilation 
and creation operators that define quasiparticle states 
will change. 

The complete set that diagonalizes the BdG Hamil- 
tonian at time t = i,(i = 1,2) is assumed to be 
{|u«),|vf),|^ ;) )}, satisfying 



i = E 



,W\/.,W 



(*)\/,,M|A. 



|uf)(uf|f 3 -|vf)(vf|f 3 



W\/,,-W| 



KW 



(64) 



As for the condensate wave function, we have normalized 
the state \<Pq ) at each time t = U as 



Vo ] (U,x):= jL= *(ti,x), 



2N, 



(0 



with 



NJf* := I n (i,:,x)d 3 a;. 



(65) 



(66) 



The particle number of the condensate will change 
through dynamic evolution. Nq is the particle number 
at each time and is a constant. 

The quasiparticle field operator can be expanded as 



/ 


« N 


\ , s 1 


M \ 1 


*,(*) 1 


u k 


1 -L A* | 


~ v k | 


Ik 


„,(*)* 


+ 7 fc 


n,W* 


\ 


-v k , 


/ \ 


, U k )\ 






(67) 



for i = 1,2. From the orthonormal relation at t = t%, we 
obtain the following expansions: 



-(2) 

Tk 



and 

~(2)t 

Tk 



+ (^ ) -^ 1)t )(ui 2) |f3|^ 1) >, 



(68) 



'7 



-E^^l^aluf 

/.(l) .(l)t\/ (2)i- I (1)\ 

-(7o -7o )( V L '\ram >• 



C,(l)t/, r (2)|- ,„(1) 



k F3|V,- 



(69) 



It can easily be shown that the bra-kets in the preceding 
expressions satisfy the following relations: 



(4 2) |T3 



(1) ) 



(u 



(2) I- I (1)\ 



,( 2 ) 



,( 2 ) 



1"3|V 



(!)\ 



73 U 



W\ 



/ (2)|- I (1)\ / (2)i- I 1 



(70) 
(71) 
(72) 



Since we have assumed that the initial condensate 
has no excitations, the quasiparticle state is the ini- 
tial vacuum denoted by |0)(!), satisfying j k |0)(!) = 0. 
Then, the expectation value of the final number operator 
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The coefficient of the second term on the right-hand side 
vanishes for the following reason. We have assumed that 
the system is quasi-static at t = tz, then the ortho normal 
relation between states implies that 



In the remainder of this section, we discuss the re- 
lation to the usual Bogoliubov theory, which does not 
deal with the backreaction to condensates. To reduce the 
usual Bogoliubov theory, we need to neglect the correla- 
tion terms {<jy<j>) and {<f>4>) in the extended GP equation 
(Eq. (|3])) and in the BdG equation (Eq. ©). The usual 
Bogoliubov theory is also gapless, and the condensate 
wave function is a solution of the BdG equation. In the 
usual Bogoliubov theory, the particle creation spectrum 
can be derived in exactly the same manner as that given 
in this section. Of course, in the Bogoliubov theory, the 
particle number of the condensate is conserved, implying 
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This point is different, but it does not af- 
fect the result. The expression for the particle creation 
spectrum given by Eq. (|77p is also valid for the usual 
Bogoliubov theory. 

The expression of the particle creation spectrum (Eq. 
([77])') agrees with that obtained in Ref. J25J, which is 
calculated using the analogy with QFT in curved space- 
time. Therefore, the particle creation discussed in Ref. 
|25| can be thought of as a phenomenon explained by re- 
diagonalization of the BdG Hamiltonian within the usual 
Bogoliubov theory. 
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Using Eqs. (f74|) and (|75|). it can be shown that 
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If the right-hand side of Eq. ((77)) is not zero, quasipar- 
ticles will appear at t = ti- This is the formula for the 
particle creation spectrum in Kita theory. The spectrum 
given by Eq. ([77]) includes the backreaction effect. This 
expression agrees with Eq. (|6Tj) calculated by QFT in 
effective spacetime. 

To evaluate the integral (|78|) at t = t%, one has to solve 
the time-dependent BdG equation to obtain the value of 

functions u k and v k at t = ti- Therefore, particle cre- 
ation reflects the history of how the condensate evolves. 



We have shown that particle creation in BEC can be 
formulated as rediagonalization of the BdG Hamiltonian 
in Kita theory. The expression for the particle creation 
spectrum has been obtained, where the quantum backre- 
action effects are taken into account. Surprisingly, it is 
the same, apparently, as one obtained within Bogoliubov 
theory [25|], but the wave functions in Eq. tfTB"]) obey the 
equation affected by the backreaction and are quantita- 
tively different from those in Bogoliubov theory J25J, in 
general. The analogy with curved spacetime has been 
constructed within Kita theory. By using QFT in effec- 
tive spacetime, we have obtained an expression for the 
particle creation spectrum, which gives the same spec- 
trum as that obtained by rediagonalizing the BdG Hamil- 
tonian. 

This implies that the usual prescription for particle cre- 
ation in QFT in curved spacetime might give the correct 
result. In the analogy, the explicit backreaction effect 
appears only in the effective spacetime metric. 

The BdG Hamiltonian rediagonalization formulation 
provides a rigorous result for particle creation from a 
theoretical point of view, in the sense that the hydrody- 
namic approximation (which is assumed in constructing 
the analogy) is not required in this formulation. Recently, 
Barcelo, Garay and Jannes j26| formulated the analogy 
between Bogoliubov quasiparticles and quanta in curved 
spacetime beyond the hydrodynamic approximation by 
introducing a generalized Klein-Gordon inner product. 
Using this inner product, the calculation based on QFT 
in curved spacetime will also give a rigorous result in the 
high energy regime. 



In the BdG rediagonalization analysis, we assumed 
that the quasiparticle vacuum (|0)(i)) is time indepen- 
dent, as is usual in QFT in curved spacetime. If the 
quasiparticle vacuum evolves in the same manner as the 
condensate and, at the final time, the initial vacuum be- 
comes the final vacuum (|0)(i) = |0)(2)), then particle 
creation will not occur in the BEC system, at least at 
zero temperature [32J . In Kita theory, it is not clear how 
to represent the quasiparticle vacuum [39J . To verify this 
assumption, a deeper understanding of the quasiparticle 
vacuum and the Bose condensate state is required, in 
addition to further experimental investigation. 

In this article, we used Kita theory as a conserving 
gaplcss mean field theory. There is another theory for- 
mulated by Yukalov which also satisfies conservation laws 
and has gapless excitations J3EH37I] . On the basis of this 
theory, particle creation with quantum backreaction will 
be able to be discussed. 
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lowing relations: 
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The Hermite conjugate of Eq. (|A2[) should be given by 
Eq. (g3|, leading to 



I (2),„ , (1). , (2),~ | (1) 



(A7) 



Therefore, it must be 

e = -l. (A8) 

Appendix B: Equations for mode functions 

To derive Eq. (|46p . we write the equations for Uj + 
v* and Uj — v*. The equation for the condensate wave 
function (Eq. ([3])) leads to 



Appendix A: quasiparticle field expansion 

We consider a field expansion with two different com- 
pete sets {u« v[}\ <p™} and {u< 2 \ vf\ ^ 2) }. The 
quasiparticle field operator can be expanded in terms of 
the two complete sets as 
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for i = 1,2. We have introduced e = ±1 in the last term 
of the expansion equation (|A1|) . which is what we seek to 
determine in this Appendix. 

The orthonormal relations for these states are given by 

Eq. (fTTj) . By using these relations, we can relate 7^. and 

.(2) 

Tj a s 



.(2) 
Ik 



EbfVfNuf) 



■7. 



(l)t/„(2) 



U fe |73|Vj- 



(l)i 



j¥0 



+ (7^+^ 1)t )(ui 2) |f3|^ 1) ), 



and 



»(2)t 
Ik = 



£ W 



,< 2 



k \-ni\Uj 



(!)\ 



^Wt/,^ 2 ) 



1 



k l T 3|V 



(A2) 



(i)\ 



i#o 



-(^ 1) +^ 1)t )(vf ) |f 3 |^). 



(A3) 



The bra-kets in the preceding expressions satisfy the fol- 
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The BdG equations (|3"5|) then give equations for the func- 



tion Uj + v* 



ih(d t + D v ){uj + v*) 
2n U -^((W) + (ftft))\(uj-v*) 
Uoifififa + Uoifrfovl, (B3) 



and for Uj — v* 

ih{d t + D v )(v,j — Vj) = K'(v,j + v*) 

~^o((#)-(W>)(«i-5;), (B4) 

where we have used the hydrodynamic approximation 

2n C/ + A"«27i t/ . (B5) 

The extended GP equation (Eq. [3]) leads to the follow- 
ing relations for any function F 

{d t + D V )^F = ^{dt + v • V)F 
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K'(^mF) = -|^V • (noVF). (B8) 

Eq. (|B6p can be applied to Eq. (|B3[) , leading to 
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From this equation, we obtain Eq. (|53[) . Substituting 
this expression into Eq. (|B4|) . we obtain 
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Then, Eq. (|B7j| leads 
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= K'(u j+ v*). (Bll) 

The right-hand side can be rewritten using Eq. (IB8I) 
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Therefore, the mode functions given in Eq. (J4"5)) satisfy 
the wave equation 



- (ft + V • v ) 7rl r(ft + v • V)/j + V • n V/j = 0, 

(B13) 

which can be rewritten in the form of Eq. (|46[) with the 
effective metric 
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